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(57) Abstract 

Timing recovery, with or without frame synchronization, in a cellular TDMA communications system is carried out by a method 
comprising estimating iiidirect variables, which embody sampling delay information, of a linear complex vector using a maximum likelihood 
criterion, and recovering the sampling delay from the estimated indirect variables. For recovering the sampling delay, an optimal procedure 
and three sub-optimal, but computationally simpler, procedures are described, along with their physical implementations. For frame 
synchronization, the timing recovery sampling delay is supplemented by a number of sample spacings determined by a maximum of signal 
amplitudes of a plurality of samples, calculated from the estimated indirect variables and the timing recovery sampling delay. 
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TIMING RECOVERY AND FRAME SYNCHRONIZATION IN A CELLULAR COMMUNICATIONS 
SYSTEM. 

This invention relates to tuning recovery and frame synchronization in 
communications systems. The invention is applicable to any communications system 
5 having known sync (synchronization) sequences or words, and is especially applicable to, 
and is described below in the context of, aTDMA (time division multiple access) cellular 
communications system compatible with EIA/TTA document IS-54-B: Cellular System 
Dual-Mode Mobile Station-Base Station Compatibility Standard (Rev. B). For 
convenience and brevity, such a system is referred to below simply as an IS-54 system. 

10 In such a system, data is communicated in time slots each comprising a sync word of 14 
symbols followed by an information sequence. 
Background of the Invention 

It is well known that it is necessary in communications systems to recover the 
timing and synchronize to the time division multiplex (TDM) frames of a received digital 

15 data communications signal, so that samples of the signal are obtained at optimum times 
for further processing to recover the communicated data. It is also well known that timing 
recovery, frame synchronization, and the necessary processing of the samples are made 
more difficult by a low signal-to-noise ratio (SNR), and that a low SNR can often be 
present in cellular communications systems. 

20 Denoting the symbol spacing of the received signal, and hence the period between 

successive samples, as T, and denoting the sampling delay, Le. the period between the 
optimal and actual sampling times of the received signal, as x, then frame synchronization 
serves to ensure that the sampling delay x is within one symbol spacing, i.e. x is within 
the interval from -T/2 to T/2, and timing recovery serves to reduce the sampling delay x to 

25 substantially zero. In practice, a sampling delay control signal can be used to adjust the 
actual sampling times or, equivalently, to control an interpolator to which the actual 
samples are supplied to obtain interpolated samples at the optimal sampling times, 
whereby frame synchronization and timing recovery is achieved. 

An object of this invention is to provide improved frame synchronization and 

30 timing recovery in a communications system. 
Summary of the Invention 

According to one aspect, this invention provides a method of determining 
sampling delay in samples of a received communications signal, comprising the steps of: 
estimating, using a maximum likelihood criterion, indirect variables of a linear complex 

35 vector which approximates the received signal samples, the indirect variables embodying 
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sampling delay information; and determining a sampling delay from the estimated indirect 
variables for timing recovery of the received communications signal 

Preferably each component of the linear complex vector comprises a linear 
combination of a predetermined pair of functions dependent upon the sampling delay, and 
5 the pair of functions conveniently comprise die functions cos(?rc) and sin(7cx) where % 
represents the sampling delay. Each indirect variable is conveniently estimated by filtering 
the received signal samples using a finite impulse response filter characteristic. 

The sampling delay can be determined from the estimated indirect variables in at 
least the following possible ways involving a trade-off between accuracy and 
10 computational complexity: 

0) by multiplying the estimated indirect variables by real elements of 
predetermined stored matrices; 

(ii) by performing a single non-linear transformation of the estimated indirect 
variables: 

IS (iii) by calculation in accordance with a fiinction(l/rc)atanof therealpartof a 

product of the complex conjugate of a first estimated indirect variable with a second 
estimated indirect variable, divided by the square of the amplitude of the first estimated 
indirect variable; 

Qv) by a sequence of steps comprising initially calculating an estimate of the 

20 sampling delay in accordance with a predetermined function of the estimated indirect 
variables, and iteratively one or more times: estimating a complex fading factor of the 
received signal samples in dependence upon the estimate of the sampling delay and the 
estimated indirect variables; and re-estimating the sampling delay in dependence upon the 
estimated complex fading factor and the estimated indirect variables. 

25 In order to provide for frame synchronization as well as timing recovery, the 

method can further comprise the steps of: adding to the determined sampling delay for 
timing recovery of the received communications signal a sampling delay comprising an 
integer number i of sampling spacings T for frame synchronization of the received 
communications signal; and determining the number i by the steps of: determining, in 

30 dependence upon the estimated indirect variables and the determined sampling delay for 
timing recovery, signal amplitudes of a plurality of received signal samples; and selecting 
the number i corresponding to a maximum amplitude sample. The step of determining 
signal amplitudes can comprise iteratively determining the signal amplitudes in dependence 
upon iterative determinations of the sampling delay for timing recovery. 

35 Another aspect of the invention provides a method of frame synchronization and 

timing recovery by determining sampling delay in samples of a received communications 
signal, comprising the steps of: estimating, using a maximum likelihood criterion, indirect 
variables of a linear complex vector which approximates the received signal samples, the 
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indirect variables embodying sampling delay information; and for each of a plurality of 
received signal samples: determining a first sampling delay from the estimated indirect 
variables, the first sampling delay being less than the sample spacing; determining, in 
dependence upon the estimated indirect variables and the respective first sampling delay, 
5 the signal amplitude of each of a plurality of received signal samples; identifying an index 
of a maximum of the determined signal amplitudes relative to a current sample; and adding 
the first sampling delay for the sample identified by the index to a product of the sampling 
spacing multiplied by the index to produce a combined sampling delay for frame 
synchronization and timing recovery. 

10 The steps of determining the first sampling delay and the signal amplitudes can be 

performed iteratively. 

The invention also provides apparatus for determining a sampling delay for 
sampling a received communications signal, comprising: a plurality of finite impulse 
response filters responsive to received signal samples for producing a plurality of indirect 

15 variables of a linear complex vector which approximates the received signal samples 
according to a maximum likelihood criterion; and a calculation unit responsive to the 
plurality of indirect variables for calculating an estimated sampling delay. The plurality of 
finite impulse response filters and die calculating unit are conveniently constituted by 
functions of at least one digital signal processor. 

20 Brief Description of the Drawings 

The invention will be further understood from the following description with 
reference to the accompanying drawings, in which: 

Fig. 1 schematically illustrates a block diagram of parts of a wireless digital 
communications receiver, 

25 Fig. 2 illustrates a block diagram of a timing recovery arrangement in accordance 

with this invention; 

Figs. 3, 4, and 5 schematically illustrate in greater detail timing recovery 
arrangements in accordance with embodiments of the invention; and 

Fig. 6 schematically illustrates a frame synchronization and timing recovery 

30 arrangement in accordance with a further embodiment of the invention. 
PetailpdPescription 

The following description initially presents, by way of example, signal and 
observation models for an IS-54 system, followed by a description of procedures which 
can be used in accordance with the method of the invention. Physical implementations of 

35 timing recovery arrangements and a frame synchronization and timing recovery 

arrangement in accordance with the invention are then described in detail with reference to 
the drawings. Although the detailed description relates specifically to IS-54 systems , it is 
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emphasized that this is by way of example and that the invention is applicable to other 
communications systems with known synchronization sequences. 
Signal and Observation Models 

An IS-54 system uses 7Z/4 -shifted DQPSK (differential quadrature phase shift 
5 keyed) signal symbols which can be described by the equations: 

s k = s k _!.w k , w k = (0 k + jfl k )/V2, k = 1, 2, 3, ... (1) 

where k is a positive integer identifying the symbol Sfc, Wfc is a complex value, 
representing binary information, with real and imaginary parts G k and ^ respectively, 
OkA e {-14 } (i.e. each of 9k and ^ is one of the set of values -1 and 1, Le. is either 
10 -lorl);andls k l = lforanyk0.e.theamplitudeofs k isl). The 14 conqrfex values 
which make up an IS-54 system sync word are the set { wi, W2, ... wm) with M = 14, 
and {so, si, S2, . . . sm) represents the set of signal symbols in a sync word, with so = 1. 

With sampling as is usual at twice the symbol rate, a discrete observation model of 
the received signal samples has the form: 

15 . yi=Ui(Sj 1 =o s ^( i ^- t - kT -K)) + ^i (2) 

where yi is a complex sample identified by the index i which is an integer from 1 to 
2M+1, T is the symbol spacing, x is the sampling delay, U| is an unknown complex 
fading factor, g(t) is the impulse response of the channel filters (the transmit and receive 
filters combined) given by: 

20 = ( sin(7rt/T) V cos(tt7tt/T) 

I «/T J[(l-(2at/T) 2 ) ; 



where a is the filter roll-off coefficient, and T|i is the sequence of complex Gaussian 
random variables with zero mean, variance 2otj, and correlation function 
2 <*ng(((i-j)T)/2) between two random variables t\i and T|j . 
Timing Recovery Method 
25 Hie observation model given by equation (2) is approximated by a linear vector 

observation model containing indirect variables which embody information as to the 
sampling delay x. These indirect variables are estimated according to the maximum- 
likelihood criterion. The sampling delay is then recovered from the estimated indirect 
variables. 

30 Linear Vector Observation Model 

Withf ] T representing the conjugate transpose of the( fiiarix ^vithin the square 
brackets, let the (2M + l>dimensional observation vector of the samples of the received 
signal be: 

Y n = [y2n+l*y2n+2— y2n+2M»y2n+2M+l] T » (3) 
35 the (M + l)-dimensional vector of the known symbols of the sync word be: 



WO 97/08867 



PCT/CA96/00429 



T 

S-=[8o.Si.—Sm-1»8m] » 

and tbe (2M + l)-dimensional observation vector of noise samples be: 

H n - [ T l2n+l» T l2n+2»— T l2n+2M» T l2n+2M+l] T ' 
Denoting tbe impulse response matrix by: 

gHO g(-*-T) - g(-x-MT) 

g(-x+T/2) g(-x-T/2) - g(-x-MT+T/2) 
G(x) = g(-x+T) g(-x) g(-x-(M-l)T) 



(4) 
(5) 



(6) 



g(-x+MT) g(-x + (M-l)T) — g(-x) 

tbe model of equation (2) can then be written in the vector form: 

Y n =G(x)SU n +H n (7) 
In order to linearize this model, each component gij(x) of the matrix G(x) is 
approximated by a linear combination, plus a constant term, of some function pair <pi(x) 
10 and 92(x), so that 

BijW - a^j q> 2 (x) +a2ij<P2(T)+a 3ij 

where aiij, a2ij, and asy are matrix coefficients for matrices Ai, A2, and A3 respectively. 
Then it is possible to make the following approximation for the matrix G(x): 

G(x)» A^x)* A 2 q>2(x)+ A 3 (8) 
15 in the interval from -T/2 to T/2 for the sampling delay x, 

A number of function pairs, for example represented by the following pairs of 
equations (9) to (13), can be used for this linearization and provide sufficient 
approximation accuracy. In the pair of equations (1 1), g(x) represents the Hilbert 
transform of the function g(x): 
20 <p 1 (x) = cos(7ix) q>2(x) = sin(7cx) (9) 

<p 1 (x) = cos(jcx/2) <p 2 (x) = sin(jn/2) ( 10 ) 

<Pl(x) = g(x) q>2W = K x ) (ID 

<Pl(x) = g(x) <P2(x) = dg(x)/dx (12) 

9l(x) = (g(x + T/4)+g(x-T/4))/2 <f>2(x) = (g(x + T/4)-g(x-T/4))/2 (13) 

25 While any of these or other possible function pairs can be used, the first function pair (9) 
provides good accuracy and the least computational complexity and is preferred for these 
reasons. Only this function pair is considered in further detail below. 

Using tfie approximation in equation (8), the observation model of equation (7) 
becomes 

30 Y n = (A, q>l(x)+ A 2 q>2(x)+ A 3 )SU n + H n (14) 

orequivalently: 

Y n = Aj S U n q>! + A 2 S U n <(> 2 + A 3 S U n + H n (15) 
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If is a variable which is a 3-dimensional complex vector constituted by the 
transpose of three indirect variables (J^ fc*, and fe^ Le. <X> n = [<|> 1>n , fc*. fc^jT with 
4*1^ = U n 9l» $2^1 = U n q>2, and fo^ = U n , and with B being a known constant matrix 
B - [AiS t A2S, A3S], then equation (15) can be written as: 
5 Y n =B0> n +H n (16) 

Estimation of Indirect Variables 

ff Vn is the known correlation matrix of the additive Gaussian noise vector , 
then from equation (16) it can be seen that the likelihood function for estimating the 
indirect variable vector <&h is given by the conditional probability p( Y n 1 4^ )0.e. the 
10 probability of Y n given the condition ^): 

expf-((Y n -BiI>jV^(Y n -B*jy^ 

pYn *n =— ^ ' . Ah J1J - (17) 

(2n) M+ V 2 det(v T ,) V2 

The suffix 1 indicates the conjugate transpose. Representing estimated values by a 
circumflex A , the maximum likelihood estimate &a for the indirect variable vector &r can 
be determined to be: 

15 *n=CY n (18) 

where C = ( B^-lfi^lBV^l is a 3 x (2M + 1) matrix that can be calculated from B and 
V|j (both of which are known) and stored in a look-up table in memory. 

The next step is to recover the sampling delay x from the estimate . 
Recovery of Sampling Delay 

Any of several procedures can be used to recover the sampling delay, the choice 
depending on the trade-off between estimation accuracy and computational complexity. 
The following describes an optimal estimation and three sub-optimal, but computationally 
simpler, estimations* 
Optimal Estimation 

25 From equation (17), a new observation equation for indirect variables can be 

derived: 

*n=*n+r n (19) 
where r n is a 3-dimensional vector of complex Gaussian random variables with known 
correlation matrix V y = ( BV^B)-!. Equation (19) can be written in the form 
30 *h=FWU n +r n (20) 

where F(t) = [ <pi(x), q>2(x), 1] T . The observation noise in equation (20) is Gaussian, so 
that die likelihood function of the observation can be obtained as: 

exp|-((* D -F(t)U„)'v T -l(*„-F( t )U„)lAl 



20 
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(23) 



and averaged over die complex variable U n to determine the following equation for the 
likelihood function with respect to v. 

exp|^V*F(T)^ 

Pi nl ^) = / \i/2 / f \ * \ (22) 

(2n)det(v T f det[F(T) V^ 1 F(t)jD(<> n ) 

where D{* n ) = exp(i4»' n V 7 - 1 6 n ). 

According to the maximum likelihood criterion, an optimal estimate of the 
sampling delay x should maximize this likelihood function. Consequently, it can be 
determined mat the optimal estimate tn of the sampling delay within the time interval from 
-T/2toT/2is: 

Making the following approximation: 

*n (Qi(9i (t)- <^(x))+ Q2 q>i(x)cp 2 (x)+ 03 9lW+Q4 92 Qs)K 

where Ql to Q5 are 3 x 3 expansion matrices having teal elements that can be calculated 
and stored in a look-up table in memory, then using the function pair (9) the optimal 
sampling delay x n can be computed as: 
z n = max^q^ cos(2tct) + sin(27rc) + q^ cos(ttc) + q 4tI) sin(Tcx) + q^ ) (25) 

where q i>n = *' D Qj 4> n for i = 1, 3, 4, and 5 and q^ = *n <h&nfe- 

Equation (25) comprises a non-linear maximization, for which reason its 
implementation in practice may be computationally complex. The following three 
sub-optimal alternative procedures avoid this disadvantage and can be easily implemented 
One-Step Non-Linear Transformation 

The estimations of die three indirect variables can be expressed in the form 
$U = U n <pi(x) + for i = 1, 2, and 3, where Afc^ is the estimation error induced by 
equation (18). If the estimation error is small, then the estimated sampling delay x n can be 
approximated by a one-step non-linear transformation of the form: 



(24) 



x n =f 



(26) 



(real(fo,afe.n); 

where the function fQ depends upon the function pair 91 and q>2- Using the function pair 
(9), equation (23) becomes: 

x n =— atan ^ — t— f 

* ^reaH^^) 



(27) 
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An implementation of this first solution for the sampling delay is described below 
with reference to Kg. 3 and requires 6M + 5 complex multiplications, 6M complex 
additions, 1 real division, and 1 non-linear scalar transformation. 
Recursive Pmcediim 

The following, alternative, recursive estimation procedure for the sampling delay 
is based on a maximization of the likelihood function with respect to variables <p i , 92, and 
U n , first considering the maximum likelihood estimation for the variables 91 and 92 with 
the assumption that U n is known, and then considering the estimationofU n with the 
assumption that the sampling delay x is known. 

First, the following likelihood function can be derived from equation (21): 



with 



Hi* 



(28) 



V <P.n ~ 



V Y 11-|V Y1 3| 2 /V Y 33 



'Y22 



-1 



0 V 
= ♦ljiUn" 1 -c Y (fcu " Unju,," 1 , and = $2. n U„ 
where Vyy are elements of the matrix V Y and Cy = Vyi3 / V^. In the derivation it is 
assumed that V Y i2 = Vy2i = Vy23 = Vy^ = 0. 

Rrom equation (28) it follows that $1 j, = m fi and 92* = MZn> ^ *e estimated 
sampling delay is given by the following non-linear transformation: 



t„ =f atari 



real 



■)JJ 



(29) 



Second, from equation (21) the likelihood function for U„ assuming that the 
sampling delay x is known can be derived as: 

p($l*.knH~^(4 (W ^ (30) 
where V u =(f(x)V y " 1 F(x))" 1 and W D =^F(x)V T - 1 F(x))" 1 F(x)V Y - 1 * n , 
from which it follows that the estimate fr n = W n . 

The recursive procedure is then summarized as comprising the following four 
sequential steps for each iteration k: 

Step 1: Initial estimation of the sampling delay in accordance with equation (31): 



r r 
a tan 



rcal-)r — r 



= f 



a tan 



real 






M 2 



J) 



(31) 
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(32) 



(33) 



Step 2: Estimation of the fading factor in accordance with equation (32): 
Step 3: Re-estimation of the sampling delay in accordance with equation (33): 

{ { K. -<=,(♦*■■ -"£))„ 

Step 4: Stop if a predetermined maximum number of iterations has been reached, 
otherwise return to Step 2. 

An implementation of this recursive procedure for the sampling delay is described 
below with reference to Kg; 5. 
Simple Procedure 

A simple variant of the above recursive procedure is to stop after Step 1, without 
proceeding to Steps 2 to 4. hi this case, the sampling delay estimation is given by: 



4t 

4>2,n 



Ci 
C2 



L n> 



X„= — 



a tan 



real 



4>2,n4>l,n ; 



J) 



(34) 



where Q is die i-th element of the known matrix C 

An implementation of this simple procedure for die sampling delay is described 
below with reference to Fig, 4. 
Frame Synchronization 

The tinting recovery methods described above apply when the sampling delay x is 
in the interval from -T/2 to T/2. Frame synchronization serves to meet this requiranent 
where the absolute or total sampling delay Xs is outside of the interval from -T/2 to T/2 but 
is within a number L of symbol spacings T. Thus x$ = x + iT, where x is in the interval 
from -T/2 to T/2 and i is an index in the range from -(L- 1)/2 to (L - l)/2 where L is odd. 
Rrame synchronization serves to check all of the intervals from -T/2 + fT to T/2 + FT and 
to make a maximum likelihood choice from these intervals, Le. to make a maximum 
likelihood selection of the index L 

Frame synchronization is based on an observation interval corresponding to the 
uncertainty of L symbol spacings of die sampling delay t, so that it is applied to a set of 
observation vectors Y each of which has the form of equation (2),the set being: 

Y l£(L-lj = { Y n-(L-l)* Y n-(L»l)+2'- Y n»- Y n+{L-1)} ( 35 > 
Applying maximum likelihood principles, the task of frame synchronization and 
timing recovery jointly is to determine the maximum of the likelihood function: 
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which can be divided into the task of frame synchronization of determining the maximum 
of the average likelihood function: 

T/2 
-T/2 

and the task of timing recovery of determining the maximum of the likelihood function: 



where i = max 



Instead of directly computing die likelihood function in equation (37), an indirect 
approach can be used as follows. For each value of the index i in die range from 
-<L - 1)/2 to (L - l)/2, equations (7), (16), (18), and (20) can be expressed as: 
Y n+2 i=G(x)SU+H n+2i Y n+2i =B* + H n+2i 

*n+2i = C Y n+2i 6 n+2i = F(t)U + r n+2i 

from which the following likelihood function can be derived: 

T/2 

P(*n+2i|i)=™ J Jp(*„ +2i |x,U,i)dx.dU (39) 
-T/2U 

, H-f(*n + 2i-F(x)u) / V y - 1 (* n+2i -F(T)u)1/2 > | 

with p(*n+2i|U n ,T,i) = ' ^ 1L± (40 ) 

(2jt)^det(v y ) v 

Instead of averaging over the variable x, equation (39) is simplified by using the 
estimated sampling delay Xn+2i obtained using any of the timing recovery methods 
described above. By averaging over the variable U n , the following then results: 



, . ex pf(^n + 2iV Y - 1 F(x Il+2i )V u F(x n+2 i) V^+JA] 

,( * n+2i|i) - ^H^) ^ 

where V u =(F(x n+2i ) / V T - 1 F(x n+2i )) 1 and D(6 n+2i ) = exp(^ 

Applying the recursive procedure described above for deriving the estimate ft n , the 



likelihood function in equation (41) can be written as: 















exp 













det(V u ^DjOa+ji) 
where U n+2i =^F(x n+2i )V y - 1 F(x n+2i )j V^O'Vy" 1 *^. 
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Ignoring dependence of V n on Xn+2i» the following results: 

i = max-^p^^ili)) m max" 1 ) U n+2i | 2 ) (43) 

Thus the frame synchronization task is reduced to estimation of a signal amplitude 
for each value of the index i, and a selection of that value of the index i corresponding to 
the maximum signal amplitude. 
Combined Frame Synchronization and Ti ming Recovery 

Combining this frame synchronization method with the recursive procedure for 
timing recovery previously described results in the following recursive procedure, which 
assumes use of the function pair (9) as above, for combined frame synchronization and 
timing recovery: 

Step 1: For each value of the index i in the range from -(L - 1)/2 to (L - l)/2, perform the 

following steps 1 A to 1C: 
1 A: Estimate the indirect variable vector in accordance with equation (44): 

4>l,(n+2i) 



*n+2i 



= CY 



n+2i 



(44) 



15 1B: 



*3,(n+2i) 

Initially (for a counter k = 1 ) estimate the signal amplitude and sampling delay in 



accordance with equations (45) and (46): 



U n+2i=<t>3 f (n+2i) 



(45) 



x n+2i - 



1 



a tan 



realfo n+2i > 



'JJ 



1G Relatively (for k = 2toa maximum number K of iterations) estimate the signal 
amplitude and sampling delay in accordance with equations (47) and (48): 



(46) 



u n+2i 



'n+2i 



ik 1 
*n+2i=- 



a tan 



^D+2i) 



(47) 



(48) 



(*l,(n+2i) ~ C Y (*3.(n+2i) ~Uj +2 i)) ; 
Step 2: For each value of the index i in die range from -(L-l)/2to(L-l)/2, determine 
the index i of the maximum estimated signal amplitude determined in Step 1, and 
determine die total sampling delay from this index and the estimated sampling 
delay determined in Step 1, in accordance with equations (49) and (SO): 

* -PU\ 2 ) (49) 



i = max 



*s = fn+2i + iT 



(50) 
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The complexity of this procedure depends on the number of iterations K and on 
the number L. The procedure is greatly simplified, with a slight reduction in accuracy, in 
a similar manna- to that for the timing recovery procedure alone by eliminatin g the 
recursion step 1C, so that Step 1 only comprises the steps 1 A and IB. The equations (44) 
to (46), (49), and (50) apply as above with K=l,and with the function pair (9) equation 
(45) becomes: 

U n +2i = $Un+2i) TO K^^n+2i) + ^n+20 sto ( 3c ^n+2i) (51) 
An implementation of this combined and simplified frame synchronization and 
timing recovery procedure is described below with reference to Fig. 6. 
Physical Implementation 

Referring now to the drawings, Fig. 1 illustrates in a block diagram parts of a 
wireless digital communications receiver, in which a wireless digital communications 
signal is supplied via an RF (radio frequency) circuit 20 of a receiver to a down converter 
22 to produce a signal which is sampled by a sampler 24, the samples being converted 
into digital form by an A-D (analog-to-digital) converter 26. The digitized samples are 
interpolated by an interpolator 28 in accordance with a recovered estimated sampling delay 
x n to produce samples Y n , at estimated optimal sampling times, for further processing. 
As an alternative to the provision of the interpolator 28, the estimated sarrq>ling delay t n 
could be used directly to control the san^ling time of the sampler 24. The interpolator 28 
forms part of digital circuits 30, conveniently implemented in a DSP (digital signal 
processor) integrated circuit, which also include a timing or clock recovery and frame 
synchronization block 32 which produces the estimated sampling delay x n as described 
below, and a carrier recovery block 34 which is not described further here. The samples 
Y n from the interpolator 28 are supplied as the input signal to the blocks 32 and 34, 

The block 32 desirably combines the functions of frame synchronization and 
timing recovery using the procedures described above and in a manner which is further 
described below with reference to Fig. 6. However, it is also possible for the timing 
recovery procedure to be implemented separately from the frame synchronization 
procedures, and to this end initially only the timing recovery functions are described 
below with reference to Figs. 2 to 5, it being assumed in these cases that the sampling 
delay tn is.within the interval from -T/2 to T/2, 
Timing Recovery 

Fig. 2 illustrates a block diagram of an implementation of the timing recovery parts 
of the block 32 for producing the estimated sampling delay x n from the received signal 
samples Y n in accordance with the optimal estimation procedure described above. This 
implementation comprises a FIR (finite impulse response) filter unit 40, a multiplier unit 
42, and a calculation unit 44. The FIR filter unit 40 is supplied with the signal samples Y n 
and filters these to produce the estimate &n in accordance with equation (18) above. The 
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multiplier unit 42 produces the product &nQi&n in accoidance with equation (24) above 
from the estimate <^ and matrices Qj, for values of i from 1 to 5, obtained as described 
above from look-up tables, the output of the multiplier unit 42 comprising constant scalars 
qj. These are supplied to the calculation unit 44 which produces the estimated sampling 
5 delay t n in accordance with equation (25). 

Fig. 3 illustrates a block circuit diagram of an implementation of the timing 
recovery parts of the block 32 for producing the estimated sampling delay x n from the 
received signal samples Y n in accordance with the one-step non-linear transformation 
procedure described above* In this implementation, the signal samples Y„ are supplied to 

10 three FIR filters 51, 52, and 53 for producing the indirect variable estimates $i, fet and $3 
respectively in accordance with equation (18). The real and imaginary parts of these 
estimates are separated by units 54, 55 , and 56 respectively and are used by the remainder 
of the circuit of Kg. 3 to implement the function of equation (27). Hiis part of the circuit 
comprises two complex multipliers 57 and 58, each comprising two multipliers for 

15 multiplying real and imaginary components and an adder for summing the products of 
these multipliers, arranged to produce respectively the denominator and the numerator in 
equation (27), a divider 59 arranged to perform fee division of the numerator by the 
denominator, and a calculation unit 60 arranged to perform the function (l/jc)atan() of 
equation (27) and hence to produce the estimated sampling delay x n . 

20 Rg. 4 illustrates a block circuit diagram of an implementation of the timing 

recovery parts of the block 32 for producing the estimated sampling delay x n from the 
received signal samples Y n in accordance wife the simple procedure described above. The 
same references are used as in Fig. 3 to denote similar parts. In Fig. 4, fee signal samples 
Y n are supplied to two FIR filters 51 and 52 for producing fee indirect variable estimates 

25 91 and $2 respectively in accoidance with equation (18). The real and imaginary parts of 
these estimates are separated by units 54 and 55 respectively and are used by the 
remainder of the circuit of Fig. 4 to implement fee function of equation (34). This part of 
the circuit comprises a complex squarer 61 (comprising two squarers and an adder) and a 
complex multiplier 62 (comprising two multipliers and an adder), arranged to produce 

30 respectively fee denominator and the numerator in equation (34), a divider 59 arranged to 
perform fee division of fee numerator by fee denominator, and a calculation unit 60 
arranged to perform fee function (l/jt)atan() of equation (34) and hence to produce fee 
estimated sampling delay x n . 

Rg. 5 illustrates a block circuit diagram of an implementation of fee timing 

35 recovery parts of the block 32 for producing fee estimated sampling delay x n from the 

received signal samples Y n in accordance with the recursive procedure described above. It 
incorporates the circuit of Rg. 4, shown as a unit 64, for producing fee indirect variable 
estimates $\ and (fe and fee estimated sampling delay Xn (Le. fee initial estimated sampling 
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delay in Step 1 of the recursive procedure, for which k = 0), and the FIR filter 53 of 
Fig. 3 arranged to produce the indirect variable estimate $3. The circuit also comprises 
units 65 to 71 and a switch 72 whose functions are described below. 

Step 1 of the recursive procedure is earned out by the unit 64 as described above 
5 to provide the initial estimated sampling delay for k = 0, thereby implementing equation 
(3 1) (or (34)), the switch 72 being in the position shown to supply this estimate to the 
output and to the calculation unit 65. The calculation unit 65 is also supplied with the 
indirect variable estimates §1 and and the matrix Vy, and calculates the estimate On 
(Step 2 of the recursive procedure) in accordance with equation (32). The units 66 to 71 

10 implement equation (33), corresponding to Step 3 of the recursive procedure, to produce 
the estimated sampling delay for the next-higher value of k, for which k > 0 so that as 
indicated in Rg. 5 the switch 72 is moved to its other position to provide the new 
estimated sampling delay to the output and to the calculation unit 65. 

As can be appreciated from Fig. 5 and equation (33), the unit 66 is an adder, with 

15 a subtraction input, which performs the subtraction in brackets in the denominator of 

equation (33), the difference is multiplied by Cy in the multiplier unit 67 and the product is 
subtracted from <|>x in the adder unit 68 to produce the denominator in equation (33). The 
unit 69 is a divider which performs the division in equation (33), die unit 70 provides the 
real part of the division result, and die unit 71 is a calculation unit which performs the 

20 function (l/jc)atan0- It can be appreciated that the calculation units 71 and 60 (in the unit 
64) can be constituted by a single unit by interchanging the positions of these units and the 
switch 72. 

Timing Recovery Simulation Remits 

Computer simulation results of the different procedures described above are 

25 indicated in the table below by way of example for an SNR of 8 dB. The choice of sync 
word, from the six sync words specified in IS-54, was found to have no significant effect 
on these results. The simulation used the function pair of equation (9) in the 
approximation of equation (8). The coefficients of the constant approximating matrices Aj 
(i = 1 to 3) of equation (8) and of the constant expansion matrices Qi (i = 1 to 5) of 

30 equation (24) were computed separately. For each of the procedures described above, the 
table indicates the standard deviation of the estimated sampling delay (a /T), the 
probability of sampling error outside of the intervals from -0.15T to 0.15T, -0.2T to 
0.2T, and -0.25T to 0.25T, denoted P.15, P.2, and P.25 respectively, and the number of 
floating point operations, denoted Nflops, for each simulation. With other SNRs, the 

35 latter number is largely unchanged 
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Optimal 
Estimation 


Non-linear 
Transformation 


Simple 
Procedure 


Recursive 
Procedure 




Kg. 2 


Fig. 3 


Fig. 4 


Fig. 5 


Delay Equation: 


(25) 


(27) 


(34) 


(33) 


a/T 


0.0426 


0.0477 


0.051 


0.0462 


P.15 


0.0008 


0.0019 


0.0035 


0.0012 


P.2 


0 


0.0001 


0.0003 


0 


P.25 


0 


0 


0 


0 


Nflops 


1461 


710 


480 


893 



As can be seen from the table, the optimal estimation provides the best results, at 
the expense of computational complexity (high Nflops). The sub-optimal procedures 
5 involve considerably less computational complexity, with the recursive procedure of 
equations (31) to (33) illustrated by Fig. 5 providing a presently preferred compromise 
between accuracy and computational complexity. 
Combined Frame Synchronization and Uminf? Recovery 

Fig. 6 illustrates a block circuit diagram of an implementation of combined frame 

10 synchronization and timing recovery functions of the block 32 for producing the estimated 
sampling delay Xg from the received signal samples Yn+2i in accordance with the first step 
of the recursive procedure described above. It incorporates the circuit of Fig. 4, again 
shown as a unit 64, which produces the real (Re) and imaginary (Im) components of the 
indirect variable estimates <|>i t (n+2i) and <|>2 f (iH-2i) as well as the estimated sampling delay 

15 x n +2i on a line 72 as described above and in accordance with equation (46). 

The circuit of Fig. 6 also comprises units 73 to 83 which serve to produce on a 
line 84 squared estimated signal amplitudes I frn+2i I 2 in accordance with equation (51). 
The units 73 and 74 comprise cos and sin operators (e.g. look-up tables in memory) 
which are supplied with the estimated sampling delay x n+ 2i from the line 72 and produce 

20 at their outputs the values cos(7tXiH.2i) and sin(7ix n+ 2i) respectively, used in equation (51). 
The units 75 to 78 are multipliers which multiply die real and imaginary components of the 
indirect variable estimates <J>i,(n+2D and fcftHZi) from the unit 64 by the outputs of the 
units 73 and 74 to implement the products in equation (51). The units 79 and 80 are 
adders arranged to add together respectively the real and imaginary components of these 

25 products, the units 81 and 82 are squarers arranged to square die outputs of the adders 79 
and 80 respectively, and the unit 83 is an adder arranged to add together the outputs of the 
squarers 8 1 and 82 thereby to produce the squared estimated signal amplitudes I tj^^i f 2 
on the line 84. 
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In the remaining parts of the circuit of Kg. 6, the estimates on the line 72 are 
sampled by a sampling switch 85 at the sampling times iT and the samples are supplied to 
an L-stage shift register 86, and similarly the estimates I frn+a I 2 on the line 84 are 
sampled by a sampling switch 87 at the sampling times iT and the samples are supplied to 
5 an L-stage shift register 88. For example, L = 5. A unit 89 is arranged to determine, at 
each sampling time and in accordance with equation (49), the index i of the maximum 
squared amplitude in the L stages of the shift register 88, and supplies this index as an 
output to a multiplier 90 and to a selection control input of a selector 91. The selector 91 
is supplied with the estimates x^a stoned in tire L stages of the shift register 86 and is 

10 arranged to supply to a line 92 that one of these estimates which corresponds to the index 
i at its selection control input The multipHa- 90 multq)lies the index i by the symbol 
spacing T, and an adder 93 is arranged to add the product iT to the output of the selector 
91, to produce the total estimated sampling delay % in accordance with equation (50). 
Although die above description relates only to the simple procedure for frame 

15 synchronization, it can be appreciated that this can be extended to the recursive procedure 
in a similar manner to that described above in relation to Fig. 5 for the timing recovery 
alone. It can also be appreciated that either of the procedures for frame synchronization 
can be combined with any of the procedures for timing recovery described above. 
Combined Frame Synchronization and Timinp Recovery Simulation Results 

20 Compute simulations similar to those described above for the timing recovery 

alone are illustrated in the table below by way of example for an SNR of 8 dB. The table 
compares die results for the combined frame synchronization and timing recovery 
procedure as described above, using the function pair (9) in the approximation of equation 
(8), with results for a conventional matched filter arrangement known in the art, indicating 

25 for zero and 300 Hz earner frequency shifts f the standard deviation of the estimated 
sampling delay (a / T) and the sampling error probabilities P.15, P.2, and R25. The 
choice of sync word was again found to have no significant effect on the results. For 
greater SNRs, the improvements of this combined frame synchronization and timing 
recover/ procedure over the conventional matched filter arrangement are even greater. 

30 



Carrier shift f: 


Combined 
Procedure 

0 


Matched Filter 
Arrangement 

0 


Combined 
Procedure 

300 


Matched Filler 
Arrangement 

300 


a/T 


0.053307 


0.151435 


0.056300 


0.150954 


P.15 


0.0066 


0.4016 


0.0074 


0.4008 


P.2 


0.0004 


0.2142 


0.0010 


0.2138 


P.25 


0.0002 


0.0694 


6 


0.0654 
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Although the invention has been described above primarily in relation to the 
function pair of equation (9), it should be appreciated that any of the other function pairs 
in equations (10) to (13), or other function pairs for use in the approximation of equation 
(8), may be used. In addition, although particular embodiments of the invention have 
5 been described in detail, it should be appreciated that numerous other modifications, 

variations, and adaptations may be made without departing from the scope of the invention 
as defined in the claims. 
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WHAT IS CLAIMED IS: 

1 . A method of determining sampling delay in samples of a received communications 
signal, comprising the steps of: 

estimating, using a maximum likelihood criterion, indirect variables of a linear 
5 complex vector which approximates the received signal samples, the indirect variables 
embodying sampling delay information; and 

determining a sampling delay from the estimated indirect variables for timing 
recovery of the received communications signal 

2 . A method as claimed in claim 1 wherein each component of die linear complex 
1 0 vector comprises a linear combination of a predetermined pair of functions dependent 

upon the sampling delay. 

3. A method as claimed in claim 2 wherein the pair of functions comprises functions 
cos(tcx) and sin(7rr) where % represents the sampling delay. 

4 . A method as claimed in any of claims 1 to 3 wherein the step of estimating each 
15 indirect variable comprises filtering the received signal samples using a finite impulse 

response filter characteristic. 

5. A method as claimed in any of claims 1 to 4 wherein the step of determining the 
sampling delay from the estimated indirect variables comprises multiplying the estimated 
indirect variables by real elements of predetermined stored matrices. 

20 6. A method as claimed in any of claims 1 to 5 wherein the stq) of determining the 
sampling delay from the estimated indirect variables comprises performing a single non- 
linear transformation of the estimated indirect variables. 

7. A method as claimed in any of claims 1 to 5 wherein the step of determining the 
sampling delay from the estimated indirect variables comprises calculating the sampling 

25 delay in accordance with a predetermined function of the realpartof a product of the 
complex conjugate of a first estimated indirect variable with a second estimated indirect 
variable, divided by die square of the amplitude of the first estimated indirect variable. 

8 . A method as claimed in claim 7 wherein the predetermined function comprises the 
function (l/7c)atanQ- 
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9. A method as claimed in any of claims 1 to 5 wherein the step of determining the 
sampling delay from the estimated indirect variables comprises the steps of initially 
calculating an estimate of the sampling delay in accordance with a predetermined function 
of the estimated indirect variables, and iteratively one or more times: 

5 estimating a complex fading factor of the received signal samples in dependence 

upon the estimate of the sampling delay and the estimated indirect variables; and 

re-estimating the sampling delay in dependence upon the estimated complex fading 
factor and the estimated indirect variables. 

10. A method as claimed in any of claims 1 to 9 and further comprising the steps of: 
10 adding to the determined sampling delay for timing recovery of the received 

communications signal a sampling delay comprising an integer number i of sampling 
spacings T for frame synchronization of the received communications signal; and 
determining the number i by the steps of: 

determining, in dependence upon the estimated indirect variables and the 
15 determined sampling delay for timing recovery, signal amplitudes of a plurality of received 
signal samples; and 

selecting the number i corresponding to a maximum amplitude sample. 

11. A method as claimed in claim 10 wherein the step of determining signal amplitudes 
comprises iteratively determining the signal amplitudes in dependence upon iterative 

20 determinations of the sampling delay for timing recovery. 

12. A method of frame synchronization and timing recovery by determining sampling 
delay in samples of a received communications signal, comprising the steps of: 

estimating, using a maximum likelihood criterion, indirect variables of a linear 
complex vector which approximates the received signal samples, the indirect variables 
25 embodying sampling delay information; and 

for each of a plurality of received signal samples: 

determining a first sampling delay from the estimated indirect variables, the first 
sampling delay being less than the sample spacing; 

determining, in dependence upon the estimated indirect variables and the respective 
30 first sampling delay, the signal amplitude of each of a plurality of received signal samples; 

identifying an index of a maximum of the determined signal amplitudes relative to 
a current sample; and 

adding the first sampling delay for the sample identified by the index to a product 
of the sampling spacing multiplied by the index to produce a combined sampling delay for 
35 frame synchronization and timing recovery. 
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13. A method as claimed in claim 12 wherein steps of determining the first sampling 
delay and the signal amplitudes arc performed iteratively. 

14. Apparatus for determining a sampling delay for sampling a received 
communications signal, comprising: 

5 a plurality of finite impulse response filters responsive to received signal samples 

for producing a plurality of indirect variables of a linear complex vector which 
approximates the received signal samples according to a maximum likelihood criterion; 
and 

a calculation unit responsive to the plurality of indirect variables for calculating an 
10 estimated sampling delay. 

15. Apparatus as claimed in claim 14 wherein* the calculation unit comprises a 
multiplier for multiplying the plurality of indirect variables by real elements of 
predetermined matrices from a store. 

16. Apparatus as claimed in claim 14 wherein the calculation unit comprises functions 
15 for calculating the estimated sampling delay as a fiuictionof rcal($2, $3>/real($i, $3), 

where $1, <j>2, and (fearethe plurality of indirect variables and $3 is the complex conjugate 
of $3. 

17. Apparatus as claimed in claim 14 wherein the calculation unit comprises functions 
for calculating the estimated sampling delay as a function of atan( real($2, $1) / 1 $1 1 2 ), 

20 where <|>i and $2 are the plurality of indirect variables and $1 is the complex conjugate of 

fc. ■ 

18. Apparatus as claimed in any of claims 14 to 17 wherein the calculation unit further 
comprises functions for calculating an estimated complex fading factor of the received 
signal samples in dependence upon the plurality of indirect variables and the estimated 

25 sampling delay and for recursively calculating the estimated sampling delay in dependence 
upon the plurality of indirect variables and the estimated complex fading factor. 

19. Apparatus as claimed in any of claims 14to 18 wherein die plurality of finite 
impulse response filters and the calculating unit are constituted by functions of at least one 
digital signal processor. 
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